Revealing the Formation of the Milky Way Nuclear Star Cluster via Chemo-Dynamical Modeling

The Milky Way nuclear star cluster (MW NSC) has been used as a template to understand the origin and evolution of galactic nuclei and the interaction of nuclear star clusters with supermassive black holes. It is the only nuclear star cluster with a supermassive black hole where we can resolve individual stars to measure their kinematics and metal abundance to reconstruct its formation history. Here, we present results of the first chemo-dynamical model of the inner 1 pc of the MW NSC using metallicity and radial velocity data from the KMOS spectrograph on the Very Large Telescope. We find evidence for two kinematically and chemically distinct components in this region. The majority of the stars belong to a previously known super-solar metallicity component with a rotation axis perpendicular to the Galactic plane. However, we identify a new kinematically distinct sub-solar metallicity component which contains about 7\% of the stars and appears to be rotating faster than the main component with a rotation axis that may be misaligned. This second component may be evidence for an infalling star cluster or remnants of a dwarf galaxy, merging with the MW NSC. These measurements show that the combination of chemical abundances with kinematics is a promising method to directly study the MW NSC's origin and evolution.Comment: 6 pages, 5 figures, accepted to ApJ Letter

Self-consistent modelling of the Milky Way's Nuclear Stellar Disc

© 2022 The Author(s) Published by Oxford University Press on behalf of Royal Astronomical Society. This is the accepted manuscript version of an article which has been published in final form at https://doi.org/10.1093/mnras/stac639The Nuclear Stellar Disc (NSD) is a flattened high-density stellar structure that dominates the gravitational field of the Milky Way at Galactocentric radius $30\lesssim R\lesssim 300$ pc. We construct axisymmetric self-consistent equilibrium dynamical models of the NSD in which the distribution function is an analytic function of the action variables. We fit the models to the normalised kinematic distributions (line-of-sight velocities + VIRAC2 proper motions) of stars in the NSD survey of Fritz et al., taking the foreground contamination due to the Galactic Bar explicitly into account using an $N$-body model. The posterior marginalised probability distributions give a total mass of $M_{\rm NSD} = 10.5^{+1.1}_{-1.0} \times10^8 \,{\rm M_\odot}$, roughly exponential radial and vertical scale-lengths of $R_{\rm disc} = 88.6^{+9.2}_{-6.9}$ pc and $H_{\rm disc}=28.4^{+5.5}_{-5.5}$ pc respectively, and a velocity dispersion $\sigma \simeq 70$ km/s that decreases with radius. We find that the assumption that the NSD is axisymmetric provides a good representation of the data. We quantify contamination from the Galactic Bar in the sample, which is substantial in most observed fields. Our models provide the full 6D (position+velocity) distribution function of the NSD, which can be used to generate predictions for future surveys. We make the models publicly available as part of the software package AGAMA.Peer reviewedFinal Accepted Versio

Early formation and recent starburst activity in the nuclear disk of the Milky Way

The nuclear disk is a dense stellar structure at the centre of the Milky Way, with a radius of ~150 pc (ref. 1). It has been a place of intense star formation in the past several tens of millions of years1-3, but its overall formation history has remained unknown2. Here, we report that the bulk of its stars formed at least 8 Gyr ago. After a long period of quiescence, a starburst event followed about 1 Gyr ago that formed roughly 5% of its mass within ~100 Myr, in what may arguably have been one of the most energetic events in the history of the Milky Way. Star formation continued subsequently on a lower level, creating a few per cent of the stellar mass in the past ~500 Myr, with an increased rate up to ~30 Myr ago. Our findings contradict the previously accepted paradigm of quasi-continuous star formation at the Galactic Centre4. The long quiescent phase agrees with the overall quiescent history of the Milky Way2,5 and suggests that our Galaxy's bar may not have existed until recently, or that gas transport through the bar was extremely inefficient during a long stretch of the Milky Way's life. Consequently, the central black hole may have acquired most of its mass already in the early days of the Milky Way

PHANGS-JWST First Results: A combined HST and JWST analysis of the nuclear star cluster in NGC 628

We combine archival HST and new JWST imaging data, covering the ultraviolet to mid-infrared regime, to morphologically analyze the nuclear star cluster (NSC) of NGC 628, a grand-design spiral galaxy. The cluster is located in a 200 pc x 400 pc cavity, lacking both dust and gas. We find roughly constant values for the effective radius (r_eff ~ 5 pc) and ellipticity ({\epsilon} ~ 0.05), while the S\'ersic index (n) and position angle (PA) drop from n ~ 3 to ~ 2 and PA ~ 130{\deg} to 90{\deg}, respectively. In the mid-infrared, r_eff ~ 12pc, {\epsilon} ~ 0.4, and n ~ 1-1.5, with the same PA ~ 90{\deg}. The NSC has a stellar mass of log10 (M_nsc / M_Sun) = 7.06 +- 0.31, as derived through B-V, confirmed when using multi-wavelength data, and in agreement with the literature value. Fitting the spectral energy distribution, excluding the mid-infrared data, yields a main stellar population's age of (8 +- 3) Gyr with a metallicity of Z = 0.012 +- 0.006. There is no indication of any significant star formation over the last few Gyr. Whether gas and dust were dynamically kept out or evacuated from the central cavity remains unclear. The best-fit suggests an excess of flux in the mid-infrared bands, with further indications that the center of the mid-infrared structure is displaced with respect to the optical center of the NSC. We discuss five potential scenarios, none of them fully explaining both the observed photometry and structure.Comment: 26 pages, 10 figures, 6 tables. Accepted for publication by ApJ

On the origin of a rotating metal-poor stellar population in the Milky Way Nuclear Cluster

We explore the origin of a population of stars recently detected in the inner parsec of the Milky Way Nuclear Cluster (NC), which exhibit sub-solar metallicity and a higher rotation compared to the dominant population. Using state-of-the-art N-body simulations, we model the infall of massive stellar systems into the Galactic center, both of Galactic and extra-galactic origin. We show that the newly discovered population can either be the remnant of a massive star cluster formed a few kpc away from the Galactic center (Galactic scenario) or be accreted from a dwarf galaxy originally located at 10-100 kpc (extragalactic scenario) and that reached the Galactic center

On the origin of a rotating metal-poor stellar population in the Milky Way Nuclear Cluster

We explore the origin of a population of stars recently detected in the inner parsec of the Milky Way Nuclear Cluster (NC), which exhibit sub-solar metallicity and a higher rotation compared to the dominant population. Using state-of-the-art N-body simulations, we model the infall of massive stellar systems into the Galactic center, both of Galactic and extra-galactic origin. We show that the newly discovered population can either be the remnant of a massive star cluster formed a few kpc away from the Galactic center (Galactic scenario) or be accreted from a dwarf galaxy originally located at 10-100 kpc (extragalactic scenario) and that reached the Galactic center

Revealing the Formation of the Milky Way Nuclear Star Cluster via Chemo-Dynamical Modeling

The Milky Way nuclear star cluster (MW NSC) has been used as a template to understand the origin and evolution of galactic nuclei and the interaction of nuclear star clusters with supermassive black holes. It is the only nuclear star cluster with a supermassive black hole where we can resolve individual stars to measure their kinematics and metal abundance to reconstruct its formation history. Here, we present results of the first chemo-dynamical model of the inner 1 pc of the MW NSC using metallicity and radial velocity data from the KMOS spectrograph on the Very Large Telescope. We found evidence for two kinematically and chemically distinct components in this region. The majority of the stars belong to a previously-known super-solar metallicity component with a rotation axis perpendicular to the Galactic plane. However, we identify a new kinematically distinct sub-solar metallicity component which contains about 7% of the stars and appears to be rotating faster than the main component with a rotation axis that may be misaligned. This second component may be evidence for an infalling star cluster or remnants of a dwarf galaxy, merging with the MW NSC. These measurements show that the combination of chemical abundances with kinematics is a promising method to directly study the MW NSC's origin and evolution

The Star Formation History of the Milky Way’s Nuclear Star Cluster

We report the first star formation history study of the Milky Ways nuclear star cluster (NSC), which includes observational constraints from a large sample of stellar metallicity measurements. These metallicity measurements were obtained from recent surveys from Gemini and the Very Large Telescope of 770 late-type stars within the central 1.5 pc. These metallicity measurements, along with photometry and spectroscopically derived temperatures, are forward modeled with a Bayesian inference approach. Including metallicity measurements improves the overall fit quality, as the low-temperature red giants that were previously difficult to constrain are now accounted for, and the best fit favors a two-component model. The dominant component contains 93% ± 3% of the mass, is metal-rich ( $\overline{[{\rm{M}}/{\rm{H}}]}\sim 0.45$ ), and has an age of ${5}_{-2}^{+3}$ Gyr, which is ∼3 Gyr younger than earlier studies with fixed (solar) metallicity; this younger age challenges coevolutionary models in which the NSC and supermassive black holes formed simultaneously at early times. The minor population component has low metallicity ( $\overline{[{\rm{M}}/{\rm{H}}]}\sim -1.1$ ) and contains ∼7% of the stellar mass. The age of the minor component is uncertain (0.1–5 Gyr old). Using the estimated parameters, we infer the following NSC stellar remnant population (with ∼18% uncertainty): 1.5 × 10 ^5 neutron stars, 2.5 × 10 ^5 stellar-mass black holes (BHs), and 2.2 × 10 ^4 BH–BH binaries. These predictions result in 2–4 times fewer neutron stars compared to earlier predictions that assume solar metallicity, introducing a possible new path to understand the so-called “missing-pulsar problem”. Finally, we present updated predictions for the BH–BH merger rates (0.01–3 Gpc ^−3 yr ^−1 )

Intermediate-mass black holes in globular clusters: Observations and simulations-Update

The study of intermediate-mass black holes (IMBHs) is a young and promising field of research. If IMBH exist, they could explain the rapid growth of supermassive black holes by acting as seeds in the early stage of galaxy formation. Formed by runaway collisions of massive stars in young and dense stellar clusters, intermediate-mass black holes could still be present in the centers of globular clusters, today. We measured the inner kinematic profiles with integral-field spectroscopy for 10 Galactic globular cluster and determined masses or upper limits of central black holes. In combination with literature data we further studied the positions of our results on known black-hole scaling relations (such as M •-σ) and found a similar but flatter correlation for IMBHs. Applying cluster evolution codes, the change in the slope could be explained with the stellar mass loss occurring in clusters in a tidal field over its life time. Furthermore, we present results from several numerical simulations on the topic of IMBHs and integral field units (IFUs). N-body simulations were used to simulate IFU data cubes. For the specific case of NGC 6388 we simulated two different IFU techniques and found that velocity dispersion measurements from individual velocities are strongly biased towards lower values due to blends of neighbouring stars and background light. In addition, we use the Astrophysical Multipurpose Software Environment (AMUSE) to combine gravitational physics, stellar evolution and hydrodynamics to simulate the accretion of stellar winds onto a black hole. We find that the S-stars need to provide very strong winds in order to explain the accretion rate in the galactic center